AD-A215  198 


FILE  COP* 


RESEARCH  AND  DEVELOPMENT  TECHNICAL  REPORT 


SLCET-TR-89-5 


POLARIZATION  MATRICES  OF  QUARTZ 

Arthur  Ballato  and  Theodore  Lukaszek 
ELECTRONICS  TECHNOLOGY  AND  DEVICES  LABORATORY 

July  1989 


DISTRIBUTION  STATEMENT 

Approved  for  public  release; 
distribution  is  unlimited. 


□  1  1 1. 


.  r.  ,r  ;  ~ 
t.  l .  >  »- 


DECO  G 1389  , 


EP 


US  ARMY 

LABORATORY  COMMAND 

FORT  MONMOUTH,  NEW  JERSEY  07703-5000 


NOTICES 


Disclaimers 


The  citation  of  trade  names  and  names  of  manufacturers  in 
this  report  is  not  to  be  construed  as  official  Government 
indorsement  or  approval  of  commercial  products  or  services 
referenced  herein. 


UNCLASSIFIED 


ecurity  classification  of  this  fa 


14.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2*  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

SLCET-TR-89-5 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMBNo  0704  0188 


3.  DISTRIBUTION /AVAILABILITY  Of  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBERtSI 


6b.  OFFICE  SYMBOL  74.  NAME  OF  MONITORING  ORGANIZATION 
( If  applicable) 


7b.  ADDRESS  {City.  State,  and  ZIP  Code) 


6 j.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 

US  Army  Laboratory  Command  (if  applicable) 

Electronics  Tech  &  Devices  Lab  SLCET-MA-A 


tc  ADDRESS  ( Oty ,  State,  and  ZIP  Code) 

Electronics  Technology  and  Devices  Laboratory 

ATTN:  SLCET-MA-A 

Fort  Monmouth,  NJ  07703-5000 


84.  NAME  OF  FUNDING /SPONSORING  |  Bb.  OfflCE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

ORGANIZATION  (if  applicable) 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


11.  title  (tndude  Security  Classification) 

POLARIZATION  MATRICES  OF  QuARTZ  (U) 


12.  PERSONAL  AUTHOR(S) 

Arthur  Ballato  and  Theodore  Lukaszek 


134.  TYPE  OF  REPORT  13b.  TIME  COVEREO 

Technical  Report  from  Jan  88  to  Jan  89 


16.  SUPPLEMENTARY  NOTATION 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

ELEMENT  NO. 

NO. 

NO. 

1L 162705 

H94 

K9 

WORK  UNIT 
ACCESSION  NO 

DA303394 


14.  DATE  OF  REPORT  (Year.  Month,  Day)  15.  PAGE  COUNT 

1989  July  26 


17. 

COS  \T|  CODES  | 

FIELD 

GROUP 

SUB-GROUP  1 

18.  SUBJECT  TERMS  (Continue  on  reverie  il  necessary  and  identify  by  b'O ck  number) 

Piezoelectric  resonators;  piezoelectric  transducers; 
quartz  crystals;  acousto-optics 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

In  analytical  treatments  of  piezoelectric-acoustic  transducers,  signal  processors,  and 
resonators,  the  electromechanical  transduction  mechanism  is  most  often  expressed  in  terms 
of  the  elements  of  the  piezoelectric  [e]  or  [d]  matrices.  Molecular  interpretations  uf 
piezoelectricity,  and  especially  electrooptical  applications,  usually  involve  polarization 
as  the  preferred  variable,  and  consequently  the  alternative  [a]  and  [b]  matrices  are  of 
interest.  The  elements  of  these  latter  sets  are  calculated  for  quartz  from  measured 
elastopiezodielectic  constants  taken  from  the  literature. 


20  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  [3  SAME  AS  RPT 


224  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Arthur  Ballato 


DO  Form  1473,  JUN  86 


I  21.  abstract  security  classification 
□  otic  users  |  Unclassified 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(201)  544-2773  SLCET-MA-A 


Previous  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  this  PAGE 

UNCLASSIFIED 


CONTENTS 


Page 

INTRODUCTION  .  1 

CONSTITUTIVE  EQUATION  SETS  .  1 

RELATIONS  AMONG  MATERIAL  CONSTANTS  .  5 

CALCULATION  SEQUENCE  .  7 

EXPLICIT  FORMULAS  FOR  POINT  GROUP  32  10 

INPUT  VALUES  FOR  QUARTZ .  16 

OUTPUT  VALUES  FOR  QUARTZ  .  16 

CONCLUSIONS .  19 

REFERENCES .  19 


TABLES 

Table  Page 

1.  Symbols,  Units,  and  Definitions  ......  2 

2 .  Relations  among  Material  Constants  .....  8 

3.  Further  Relations  among  Material  Constants  ...  9 

4.  Elastopiezodielectric  Matrices  for  Point  Group  32: 

The  [e] ,  [h],  and  [a]  Sets . 11 

5.  Elastopiezodielectric  Matrices  for  Point  Group  32: 

The  [d],  [g]  ,  and  [b]  Sets . ll 

6.  Isagric  Elastic  Stiffnesses  .......  16 

7.  Piezoelectric  Stress  Constants  ......  16 

8.  Dielectric  Permittivities  at  Constant  Strain  ...  16 


?  Elastic  Stiffnesses  .........  17 

10.  Elastic  Compliances  .  17 


i  i  i 


Table 


Page 


11.  Piezoelectric  [e] ,  [h] ,  and  [a]  Values  ....  17 

12.  Piezoelectric  [d] ,  [9]/  an d  [b]  Values  ....  18 


13.  Dielectric  (eps)  Values . .  .  .  18 

14.  Dielectric  (chi)  Values . .  18 

15.  Dielectric  (bet)  Values  .  ...  18 

16.  Dielectric  (zet)  Values  .  19 


INTRODUCTION 


Electromechanical  transduction  taking  place  via  the  piezoelectric 
effect  is  characterized  phenomenologically  by  constitutive  equa¬ 
tions  that  relate  the  elastic  and  electric  variables.  These  equa¬ 
tions  take  a  variety  of  forms,  depending  upon  the  choice  of  inde¬ 
pendent  and  dependent  variables;  the  choice  is  normally  dictated  by 
the  application.  For  example,  piezoelectric  resonators  in  the  form 
of  thickness  mode  plates  are  most  easily  treated  using  the  isagric 
elastic  stiffnesses  [cE],  the  piezoelectric  stress  constants  [e], 
and  the  dielectric  permittivities  at  constant  strain  [(eps)S]. 

Various  measurement  techniques  yield  values  for  the  elements  of  a 
particular  coefficient  set  more  directly  than  those  of  another. 
The  coefficients  appearing  in  the  different  equation  sets  are, 
however,  interrelated,  so  that  once  any  one  complete  set  is  avail¬ 
able,  all  the  other  sets  of  elements  may  be  found.  The  most  accu¬ 
rate  and  precise  experimental  results  to  date  have  been  from  plate 
resonator  (resonance)  and  pulse-echo  (transit-time)  measurements. 
From  the  [cE],  [e],  and  [(eps)S)  matrices  determined  therefrom, 
those  matrices  representing  material  properties  expressed  in  the 
other  alternative  forms  may  be  calculated. 

Electroopt ical  applications  are  becoming  increasingly  important. 
So  also  are  treatments  of  piezoelectric  and  ferroelectric  phenomena 
from  the  standpoint  of  molecular  interactions.  In  both  of  these 
cases  the  constitutive  equations  using  polarization  as  the  inde¬ 
pendent  electrical  variable,  rather  than  either  electric  intensity 
or  displacement,  assume  greater  importance  than  the  sets  tradition¬ 
ally  used  for  transducer,  signal  processing,  and  resonator  applica¬ 
tions. 

In  this  report  we  give  the  complete  sets  of  linear  constitutive 
equations  relating  elastic  and  electric  fields.  For  each  equation 
set  the  numerical  values  are  computed  for  quartz,  from  the  measured 
[cE],  [e],  and  [(eps)S]  values  of  Bechmann  (Ref.  1).  Coupling  to 
the  thermal  field  is  neglected.  Rationalized  mks  units  are  used 
throughout . 


CONSTITUTIVE  EQUATION  SETS 

Symbols  and  units  for  the  quantities  employed  are  given  in  Table  1. 
In  terms  of  these,  six  constitutive  equation  sets  are  used.  Of 
these,  electric  intensity,  dielectric  displacement,  and  polariza¬ 
tion  each  appear  in  two  sets  as  an  independent  variable.  The  sets 
are,  in  compressed  matrix  notation,  as  follows.  A  prime  denotes 
transpose;  [I]  is  the  unit  matrix. 

I.  The  Piezoelectric  Stress  Constant  Set 

f T)  =  [CE]  [S]  -  [e]'  (E)  (  1) 

[15]  -  [  e]  [S]  +  [(eps)S]  (E]  (  2) 
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TABLE  1.  SYMBOLS,  UNITS,  AND  DEFINITIONS. 


QUANTITY 

UNIT 

SYMBOL/DEFINITION 

Elastic  stress 

N/m2 

[T] 

Elastic  strain 

fS] 

Electric  intensity 

V/m 

[El 

Dielectric  displacement 

C/m2 

[D] 

Dielectric  polarization 

C/m2 

[P] 

Elastic  compliance  at 
constant  [E],  [D],  [P] 

m2 /N 

[sE],  [sD],  [sP] 

Elastic  stiffness  at 
constant  [E],  [D],  [P] 

N/m2 

[cEJ,  [CD],  [cP] 

Dielectric  permittivity 
at  constant  [T] ,  [S] 

F/m 

[  (eps) T] ,  [(eps)S) 

Dielectric  constant, 
relative,  at 
constant  [T] ,  [S] 

[(  Kr)T],  [(  Kr) S ] 

=[ (eps) T]/ (eps) o, 

[ (eps) S]/ (eps)o 

Dielectric  impermeability 
at  constant  [T] ,  [S] 

m/F 

( (bet) T],  Tibet) S) 

-  [  (eps)  T]  -f-  , 

[  (eps) S ]  (  1 > 

Dielectric  impermeability, 
relative,  at  constant 
[T],  [S] 

[ (betr) T] ,  [(betr)S] 

=L (bet) T] * (eps) o , 

[ (bet) S]* (eps) o 
=  [(  Kr)T]  ’j  , 

[(  Kr) S ]  i-1' 

Dielectric  susceptibility 
at  constant  [T],  [S] 

F/m 

[ (chi) T] ,  [(chi)S] 

=[ (  Kr) T-I ] * (eps) o , 

[ (  Kr) S-I ] * (eps) o 

Dielectric  susceptibility, 
relative,  at  constant 
[T],  [S] 

[ (chir) T] ,  [ (chir) S ] 

=  [  (chi)T]/ (eps)o, 

[ (chi) S]/ (eps) o 

Reciprocal  dielectric 
susceptibility  at 
constant  [T] ,  [ S ] 

m/F 

t (zet)T] ,  [ ( zet ) S ] 

=((chi)T]  “J  , 

[  (chi ) S ] 

Reciprocal  dielectric 
susceptibility,  relative, 
at  constant  [T] ,  [S] 

[  ( zetr) T] ,  [(zetr)S] 

=  [  ( zet) T) * (eps) o , 

[ (zet) S] * (eps) o 

[e] 


Piezoelectric  stress 
constant 


C/m2 
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TABLE  1.  SYMBOLS,  UNITS,  AND  DEFINITIONS,  (continued) 


QUANTITY 

UNIT 

SYMBOL/DEFINITION 

Piezoelectric 

coefficient 

strain 

m/V  == 

C/N 

[d] 

Piezoelectric 

modulus 

stress 

N/C  = 

V/m 

[h] 

Piezoelectric 

constant 

strain 

m2/C 

[g] 

Piezoelectric 

modulus 

polarization 

V/m  = 

N/C 

[a] 

Piezoelectric 

constant 

polarizat ion 

m2/C 

[b] 

Note:  Square  brackets,  sic :  [  ],  denote  matrices. 
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II.  The  Piezoelectric  Strain  Coefficient  Set 


[S]  =  [SE]  [T]  +  [d]'  [E]  (  3) 

[D]  =  [  d]  [T]  +  [(eps)T]  [E]  (  4) 

III.  The  Piezoelectric  Stress  Modulus  Set 

[T]  =  [CD]  [S]  -  [h] «  [D]  (  5) 

[E]  =  -[  h]  [S]  +  [ (bet) S ]  [D]  (  6) 

IV.  The  Piezoelectric  Strain  Constant  Set 

rs]  =  [SDJ  [T]  +  [g] '  [D]  (  7) 

[E]  =  -[  g]  [T]  +  [ ( bet) T ]  [D]  (  8) 

V.  The  Piezoelectric  Polarization  Modulus  Set 

[T]  =  [CP]  [S]  -  [a] *  [P]  (  9) 

[E]  =  -[  a]  [S]  +  [ ( zet ) S ]  [P]  (10) 

VI.  The  Piezoelectric  Polarization  Constant  Set 

[S]  =  [sP]  [T]  +  [b]'  [P]  (11) 

[E]  =  -[  b]  [T]  +  [ ( zet) T]  [P]  (12) 


The  electric  variables  are  connected  by  the  relation 

[D]  =  (eps)o  *  [E]  ^  [P]  (13) 

where  (eps)o  is  the  permittivity  of  free  space,  defined  by 

(eps)o  *  (mu)o  *  (c)  *  (c)  =  1  ;  (14) 

(mu)o  is  the  permeability  of  free  space,  equal,  by  definition,  to  4 
PI  *  lO^-7),  and  (c)  is  the  velocity  of  light  in  vacuo  and,  also  by 
definition,  is  equal  exactly  to  2.99792458  x  108  m/s. 

From  (13)  the  expressions  for  the  remaining  electric  variables 
associated,  respectively,  with  the  six  equation  sets  (  1)  to  (12) 
may  be  found: 


m 

=  [  e]  [S] 

+  [ ( chi ) S ] 

[E] 

(15) 

fp] 

=  [  cl]  [T] 

+  [ (chi ) T ] 

[E] 

(16) 

tp] 

=  (eps)o  * 

[  h]  [S]  + 

[I  - 

(eps)o  * 

(bet) S ] 

[D] 

(17) 

[pj 

=  (eps)o  * 

[  9]  [T]  + 

[I  - 

(eps)o  * 

(bet) T] 

[P] 

(18) 
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[D] 

=  -(eps)o  * 

[  a] 

[S]  + 

[I  + 

(eps)o  * 

(zet) S] 

[P] 

(19) 

[D] 

=  -(eps)o  * 

[  b] 

[T]  + 

[I  + 

(eps)o  * 

( zet ) T ] 

[P] 

(20) 

RELATIONS  AMONG  MATERIAL  CONSTANTS 

The  material  constants  are  interrelated  by  the  following  formulas: 

[CX]  [sX]  =  [ ( eps) Y ]  [ ( bet) Y ]  =  [I]  (21) 

[ (chi) Y ]  [ ( zet ) Y ]  =  [(  Kr ) Y  -  (chir)Y]  =  [I]  (22) 

In  (21)  and  (22),  X  =  E,  D,  or  P  and  Y  -  T  or  S. 

[cD]  -  [cE]  =  [h] '  [el  =  [e] '  [h] 

=  [h] •  [ ( eps ) S ]  [h]  =  [e]'  [ ( bet ) S ]  [e] 

=  [a]'  [e  -  h  *  (eps)o]  =  [e  -  h  *  (eps)o]'  [a]  (23) 

[cP]  -  [cD]  -  [ h ] *  [a]  *  (eps)o  =  [a]'  [h]  *  (eps)o 

=  [h] *  [ ( eps ) S ]  [ ( zet ) S ]  [h]  *  (eps)o 

=  [a]'  [ ( bet ) S ]  [ (chi ) S ]  [a]  *  (eps)o 

=  [a  -  h) '  [e]  =  [e] '  [a  -  h)  (24) 

[cP]  -  [cE]  =  [a] •  [e]  =  [e] '  [a] 

=  [a]'  [ (chi ) S ]  [a]  =  [e] '  [(zet)S]  [e] 

=  [h] '  [e  +  a  *  (eps)o]  =  [e  +  a  *  (eps)o]'  [h]  (25) 

[sE]  -  [sDJ  =  [d] '  [g]  =  [g] '  [J] 

=  [d  ]  •  [ (bet) t ]  [dj  =  [g] •  [(eps)t]  [g] 

)  ’  [d  -  g  *  (eps)ol  =  [d  -  g  *  (eps)o]  1  [b]  (2G) 

[sD]  -  [sP]  -  [ b ] *  [g]  *  (eps) o  =  [ g ] '  [b]  *  (eps)o 

=  [ g ] ’  [(eps)T)  [ ( zet) T]  [g]  *  (eps)o 

=  [ b ] *  [ (bet) T]  [ (chi) T]  [b]  *  (eps)o 

=  [b  -  g] •  [d]  =  [d] •  [b  -  g]  (27) 

[sE]  -  [sP]  =  [b] •  [d]  -  [d] 1  [b] 

=  [b] •  [ (chi)T]  [b]  =  [d] •  [ ( zet ) T ]  [d] 

=  [g] 1  [d  +  b  *  (eps)o]  =  [d  +  b  *  (eps)o)'  [g]  (28) 


[ ( zet )  S ]  -  [ (zet)T]  =  [b]  [a]'  =  [a]  [b]  • 

-  [b]  [cP]  [b] •  =  [a]  [sP]  [a]'  (29) 

[ (chi)T]  -  [ (chi) S]  =  [ ( eps ) T ]  -  [(eps)S] 

=  [e]  [d] •  =  [d]  [e] ' 

=  [d]  [cF]  [d] '  =  [e]  [sE]  [e]  '  (30) 

[ (bet) S ]  -  [ (bet) T ]  =  [h]  [g] '  =  [g]  [h]' 


=  [g] 

[CD] 

[g]'  =  [h] 

[sD] 

[h]  ' 

(31) 

(e]  =  [d] 

[cE] 

= 

[ (eps) S ] 

[h]  = 

[ (chi) S  ] 

[a] 

(32) 

[d]  =  [e] 

[SE] 

= 

[ (eps) T] 

[g]  = 

[ (chi) T] 

[b] 

(33) 

[h]  =  [g] 

[CD] 

= 

[ (bet) S ] 

[e]  = 

[ (chi ) S ] 

[ (bet) S ] 

[a] 

= 

[I  -  (bet) S  * 

(eps) o]  [a] 

(34) 

[g]  =  [h] 

[sD] 

= 

[ (bet) T] 

[d]  = 

[ (chi) T] 

( (bet ) T ] 

[b] 

= 

[I  -  (bet) T  * 

(eps) o]  [b] 

(35) 

[a]  =  [b] 

[cP] 

= 

[ ( zet) S ] 

[e]  = 

[ (eps) S ] 

[ (zet)S] 

[h] 

= 

[I  +  (zet)S  * 

(eps) o]  [h] 

(36) 

[b]  =  [a] 

[SP] 

= 

[ (zet) T] 

[d]  = 

[ (eps'  T] 

[ ( zet) T ] 

[g] 

=  [I  +  (zet) T  *  (eps) o]  [g]  (37) 

Some  alternative  relations  are  the  following: 


(a  -  h]  =  [(zet)S]  [h]  *  (eps)o 

=  [ (bet) S  ]  [a]  *  (eps) o  (38) 

ib  -  g]  =  [ ( zet) T]  [g]  *  ( eps) o 

=  [(bet)T]  [b]  *  (eps)o  (39) 

[e  +  a  *  (eps)o]  =  [(eps)S]  [a]  (40) 

[d  +  b  *  (eps)o]  =  [(eps)T]  [b]  (41) 

[e  -  h  *  (eps)o]  =  [(chi)S]  [h]  (42) 

[d  -  g  *  (eps)o]  =  [(chi)T]  [g]  (43) 


Equations  (21)  to  (43)  result  from  equating  1 1 k  j  dependent  variables 
in  pairs  selected  from  equations  {  1)  to  (1?)  and  (15)  to  (20).  bach 
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pair  yields  one  equation  in  three  variables,  one  mechanical  and  two 
electrical,  or  vice  versa.  Two  other  equations  exist,  again  from 
(1)  to  (12)  and  (15)  to  (20),  that  contain  the  same  three  variables 
found  in  each  paired  equation.  One  of  these  auxiliary  equations  is 
used  to  eliminate  one  of  the  two  variables  of  the  same  kind;  the 
result  is  one  equation  in  two  variables,  one  electrical  and  one 
mechanical.  These  are  now  independent  variables,  so  the  coeffi¬ 
cients  must  vanish;  two  relations  between  the  material  coefficients 
result.  As  an  example,  (  3)  and  (  7)  both  have  [S]  as  dependent 
variable.  Equating  them  produces  one  relation  in  [T],  [E],  and 
[D];  one  of  t  ,ie  electrical  variables  must  be  eliminated.  This  is 
done  by  using  either  (  * )  or  (  8) ;  ^ach  contains  the  same  three 
variables.  If  (  8)  is  used  to  eliminate  [E],  one  obtains  [sE  -  d' 
g  -  sD]  [T]  =  [d*  (bet)T  -  g']  [D].  Therefore,  [sE]  -  [sD]  =  [d]  ' 
[g]  and  [g]  =  ((bet)T]  [d].  Use  of  (  4)  instead  of  (  8)  leads  to 
the  equations  [sE]  -  [sD]  =  [g] 1  [d]  and  [d]  =  [(cps)Tj  [g] .  There 
are  36  pairs,  six  each  equating  [S]  and  [T] ,  and  eight  each  equat¬ 
ing  [E],  [D],  and  [P].  The  72  relations  contain  many  redundancies. 
Relations  between  the  elastic,  piezoelectric,  and  dielectric  con¬ 
stants  are  shown  schemac ica l ly  in  Tables  2  and  3. 


CALCULATION  SEQUENCE 

Using  as  input  [cE],  [e],  and  [(eps)S],  one  may  compute  the  remain¬ 
ing  quantities  in  a  variety  of  ways.  The  following  sequence  i* 
typical : 


[sE]  =  [cE] 

(-1) 

[ ( bet ) S ]  - 

[ (eps)S] 

( 

-i) 

[d]  =  [e]  [ 

sE] 

[h]  =  [ (bet) S ]  [e] 

[ (eps)TJ  - 

[ (eps) S] 

= 

[e] 

[d]  1 

[ ( eps) T ]  = 

[ (eps) S 1 

+ 

[e] 

[d]  • 

[ ( bet ) T ]  = 

[ (eps) T] 

(- 

1) 

[cD]  -  [cE] 

=  [e]  1 

[h] 

[CD]  =  [cE] 

+  [e]  ' 

[h] 

[g]  =  [ (bet) T ]  [d] 

[sE]  -  [sD] 

=  [d]  1 

[g] 

[sD]  -  [sE] 

-  [d]  ' 

[g] 

[ ( bet r ) S ]  = 

[ (bet) S 

]  * 

(eps)  o 

[ (zetr)S]  = 

[ (betr)S] 

[I  ~ 

( bet 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

S]  (_1)  (57) 
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TABLE  2 .  RELATIONS  AMONG  MATERIAL  CONSTANTS. 


X 

< 
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TABLf  J.  FURTHER  RELATIONS  AMONG  MATERIAL  CONSTANTS. 


(58) 


[  ( zet) S]  -  ((zetr)]  /  (eps)o 


[(betr)TJ  =  [(bet)T]  *  (eps)o 

(59) 

[  ( zetr) T]  =  [ (betrT]  [I  -  (betr)T] 

(60) 

[  ( zet) T  =  [  ( zetr) T] 

/  (eps)o 

(61) 

[  (chi) S ]  =  [  ( zet) S ] 

(-D 

(62) 

[  (chi) T]  =  [ ( zet) T] 

(-D 

(63) 

[a]  =  [ ( zet) S ]  [e] 

(64) 

[b]  =  [ (zet) T]  [d] 

(65) 

[cP]  -  [cE]  =  [e]' 

[a] 

(66) 

[cP]  =  [cE]  +  [e] ' 

[a] 

(67) 

[cP]  -  [CD]  =  [a] 1 

[h]  *  (eps)o 

(68) 

[sE]  -  [sP]  =  [d] ' 

[b] 

(69) 

[sP]  =  [sE]  -  [d] ' 

[b] 

(70) 

[sD]  -  [SP]  =  [g] • 

[b]  *  (eps)o 

(71) 

[ (bet) S]  -  [ (bet) T] 

=  (h]  [g] ' 

(72) 

[ (chi)T]  -  [ (chi) S ] 

=  [ ( eps ) T ]  -  [ ( eps ) S ] 

(73) 

[ ( zet) S ]  -  [ ( zet ) T ] 

=  [a]  [b] • 

(74) 

A  number  of  these  relations  are  used  as  checks.  For  example, 
[(bet)S]  and  [(bet)T]  are  known  from  (45)  and  (50),  but  the  differ¬ 
ence  is  recomputed  in  (72). 

EXPLICIT  FORMULAS  FOR  POINT  GROUP  32 
Elastic : 

The  6x6  elastic  constant  portions  of  Tables  4  and  5  partition  into 
4x4  and  2x2  submatrices.  The  4x4  elastic  stiffness  and  compliance 
submatrices  are  interrelated  by  formulas  (75)  to  (93),  taken  from 
Cady  (Ref.  2): 


A  = 

S3  3 

*  (sll  +  sl2 )  - 

2 

* 

s  1 3 

*  sl3 

(75) 

B  = 

s44 

*  (sll  -  sl2 )  - 

2 

* 

si  4 

*  sl4 

(76) 

2  * 

ell 

=  s3  3  /  A  +  s4 4 

/ 

B 

(77) 

2  * 

C12 

=  S3 3  /  A  -  s4 4 

/ 

B 

(78) 
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TABLE  4.  ELASTOP I EZODI ELECTRIC  MATRICES  FOR  POINT  GROUP  32: 


THE 

(e]f 

[h]  , 

AND 

[a] 

SETS. 

11 

12 

13 

14 

00 

00 

] 

11 

00 

00 

] 

cE 

3  e* 

12 

11 

13 

-14 

00 

00 

] 

-11 

00 

00 

— 

3-  - 

] 

e 

3 (eps) S 

13 

13 

33 

00 

00 

00 

] 

00 

00 

00 

14 

-14 

00 

44 

00 

00 

J 

] 

14 

00 

00 

3 

cD 

3  h' 

00 

00 

00 

00 

44 

14 

3 

00 

-14 

00 

— 

3 - 

3 

h 

] (bet) S 

00 

00 

00 

00 

14 

66 

3 

00 

-11 

00 

11 

-11 

00 

14 

00 

00 

J 

3 

11 

00 

00 

3 

CP 

3  a' 

00 

00 

00 

00 

-14 

-11 

3 

00 

11 

00 

— 

3 - 

3 

a 

] ( zet) S 

00 

00 

00 

00 

00 

00 

3 

00 

00 

33 

66 

=  (11  - 

12)  / 

2 

=  ==== 

= - : 

=  =  = 

==  =  = 

Matrix  entries  show  only  subscripts. 


TABLE  5.  ELASTOPI EZODI ELECTRIC  MATRICES  FOR  POINT  GROUP  32: 


THE 

[d3, 

[g] , 

AND  [b]  SETS. 

11 

12 

13 

14 

00 

00 

iH 

i — 

00 

00 

3 

sE 

3  d* 

12 

11 

13 

-14 

00 

00 

]  -li 

00 

00 

— 

3 

3 

d 

3  (eps)T 

13 

13 

33 

00 

00 

00 

]  00 

00 

00 

14 

-14 

00 

44 

00 

00 

J 

3  14 

00 

00 

3 

sD 

3  g' 

00 

00 

00 

00 

44 

14*2 

]  00 

-14 

00 

— 

3 - 

3 

g 

] (bet) T 

00 

00 

00 

00 

14*2 

66 

]  00 

-11*2 

00 

11 

-11 

00 

14 

00 

00 

J  ■ 

3  ii 

00 

00 

3 

SP 

1  b' 

00 

00 

00 

00 

-14 

-11*2 

]  00 

11 

00 

— 

3 - 

3 

b 

]  (zet)T 

00 

00 

00 

00 

00 

00 

]  00 

00 

33 

66  =  (11  -  12)  *  2 


Matrix  elements  show  only  subscripts. 


cl3  =  -  sl3  /  A  ;  cl4  =  -  sl4  /  B  (79a),  (79b)  (79) 

c3 3  =  (sll  +  sl2 )  /  A  (80) 

c44  =  (Sll  -  Sl2 )  /  B  (81) 

c66  =  (ell  -  cl2)  /  2  =  S44  /  (2*2)  (82) 

K  =  c33  *  (ell  +  cl2 )  -  2  *  cl3  *  Cl3  (83) 

L  =  c44  *  (ell  -  cl2 )  -  2  *  cl4  *  cl4  (84) 

2  *  Sll  =  C3 3  /  K  +  c44  /  L  (85) 

2  *  Sl2  =  C33  /  K  -  c44  /  L  (86) 

sl3  =  -  cl3  /  K  ;  Sl4  =  -  cl4  /  L  (87a),  (87b)  (87) 

s3 3  =  (ell  +  cl2 )  /  K  (88) 

s44  =  (ell  -  cl2 )  /  L  (89) 

s66  =  (Sll  -  Sl2)  =  2  *  C4 4  /  L  (90) 

det  (4X4)  [s]  =  A  *  B  (91) 

det  (4X4)  [c]  =  K  *  L  (92) 

A*K  =  B*L  =  A*B*K*L=1  (93) 

Formulas  (75)  to  (93)  hold  for  each  set  of  constant  electrical 
conditions:  either  E,  D,  or  P  constant. 

[CD]  -  [cE]  =  [del  cDE ]  =  [ e ] *  [h]  =  [h] '  [e]  (23) 

del  cDEll  =  +  ell  hll  (94) 

del  cDE12  =  -  ell  hll  (95) 

del  CDE13  =  0  (96) 

del  cDE14  =  +  ell  hl4  =  +  hll  el4  (97) 

del  CDE33  =  0  (98) 

del  CDE44  =  +  el4  hl4  (99) 

del  CDE66  =  +  ell  hll  (100) 

[cP]  -  [cD]  =  [del  cPD]  =  [a]'  [h]  *  (eps)o 

=  [h]  '  [a]  *  (eps) o  (24) 

del  cPDll  =  (  +  all  hll  )  *  (eps)o  (101) 
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del  cPD12 


(  -  all  hll  )  *  (eps)o 


(102) 


del  CPD13  =  0  (103) 

del  cPDl 4  =  (  +  all  hl4  )  *  (eps)o 

=  (  +  hll  al4  )  *  (eps)o  (104) 

del  CPD33  =  0  (105) 

del  CPD44  =  (  +  al4  hl4  )  *  (eps)o  (106) 

del  CPD66  =  (  +  all  hll  )  *  (eps)o  (107) 

[cP]  -  [cE]  =  [del  cPE ]  =  [e]'  [a]  =  [a]'  [e]  (25) 

del  cPEll  =  +  ell  all  (108) 

del  CPE12  =  -  ell  all  (109) 

del  cPEl 3  =  0  (110) 

del  CPE14  =  +  ell  al4  =  +  all  el4  (111) 

del  CPE33  =  0  (112) 

del  CPE44  =  +  el4  al4  (113) 

del  CPE66  =  +  ell  all  (114) 

From  the  del  cl4  entries  we  have  the  ratios 


ell  /  el4  =  hll  /  hl4  =  all  /  al4. 


(115),  (115) 


[sE]  -  [sD]  =  [del  sED]  =  [d]'  [g]  -  [g]'  [d]  (26) 

del  sEDll  =  +  dll  gll  (117) 

del  sEDl 2  =  -  dll  gll  (118) 

del  SED13  =  0  (119) 

del  SED14  =  +  dll  gl4  =  +  gll  dl4  (120) 

del  SED33  =  0  (121) 

del  SED44  =  +  dl4  gl4  (122) 

del  SED66  =  +  dll  gll  *  4  (123) 


[cD]  -  [sP]  =  [g] '  [b]  *  (eps)o 


= 

[b]*  Cg]  * 

(eps) o 

(27) 

del 

sDPll  = 

( 

+  gll  bll 

)  * 

(eps) o 

(124) 

del 

SDP12  = 

( 

-  gll  bll 

)  * 

(eps) o 

(125) 

del 

SDP13  = 

0 

(12€) 

del 

SDP14  - 

( 

+  gll  bl4 

)  * 

( eps) o 

= 

( 

+  bll  gl4 

)  * 

(eps) o 

(127) 

del 

SDP33  = 

0 

(128) 

del 

SDP44  = 

( 

-i-  gl4  bl4 

)  * 

( eps) o 

(123) 

del 

SDP66  = 

( 

+  gll  bll 

)  * 

4  *  (eps)o 

(130) 

[SE]  ■ 

-  [SP]  = 

f del  sEP]  = 

[b]  1 

[d]  = 

[d] '  [b] 

(28) 

del 

sEPll  = 

+ 

dll  bll 

(131) 

del 

SEP12  - 

- 

dll  bll 

(132) 

del 

sEPl 3  = 

0 

(133) 

del 

SEP14  = 

+ 

dll  bl4  = 

+  bll  dl4 

(134) 

del 

SEP33  = 

0 

(135) 

del 

SEP44  = 

+ 

dl4  bl4 

(136) 

del 

SEP66  - 

+ 

dll  bll  * 

4 

(137) 

From  ' 

the  del 

Sl4 

entries  ' 

we  have  the 

ratios 

dll  / 

dl4 

=  gii  /  ' 

gl4  = 

bll  / 

bl4 . 

(138),  (139) 

Piezoelectric 

[d]  = 

[e] 

[SE] 

(33) 

dl4  - 

+  el4  sE44  + 

ell 

SE14  * 

2 

(140) 

dll  = 

+  ell  (sEl 1 

-  sE12 )  +  el4  sE14 

(141) 

[h]  = 

f (bet) S]  [e] 

(34) 

hl4  = 

(bet)Sll  el4 

(142) 

hll  = 

(bet)Sll  ell 

(143) 
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[g]  =  [ (bet) T]  [d]  (35) 

gl4  =  (bet) Til  dl4  (144) 

gll  =  (bet) Til  dll  ( ! 4 5 ) 

[a]  =  [ ( zet ) S ]  [e]  (36) 

al4  =  (zet)Sll  el4  (146) 

all  =  (zet)Sll  ell  ^147) 

[b]  =  [ (zet) T]  [d]  (37) 

bl4  =  (zet)Tll  dl4  (148) 

bll  =  ( zet) Til  dll  (149) 

Dielectric: 

[ (bet) Y  ]  =  [ (eps) Y]  (-1>  (2i) 

(bet)Yll  =  1  /  (eps)Yll  (150) 

(bet) Y33  =  1  /  (eps) Y33  (151) 

[(zetr)Y]  =  [  (betr)  Y]  [I  -  (betrJY]^1)  (152) 

(zet)Yll  =  1  /  ((eps)Yll  -  (eps)o)  (153) 

(zet) Y33  =  1  /  ( (eps) Y33  -  (eps)o)  (154) 

( (eps) T  -  (eps) S ]  =  [del  (eps)]  =  [e]  [d] '  = 

[ (chi)T  -  (chi) S ]  =  [del  (chi)]  =  [d]  [e] •  (30) 

del  ( eps) 11  =  del  (chi)ll  =  +  el4  dl4  +  ell  dll  *  2  (155) 

del  (eps) 33  =  del  (chi)33  =  0  (156) 

[ (bet) S  -  (bet) T]  =  [h]  [g] '  =  [g]  [h] '  (31) 

del  (bet) 11  =  +  hl4  gl4  +  hll  gll  *  2  (157) 

del  (bet) 33  =  0  (158) 

[ ( zet) S  -  ( zet) T]  =  [del  (zet)]  =  [a]  [b] ’  =  [b]  [a]'  (159) 

del  (zet)ll  =  +  al4  bl4  +  all  bll  *  2  (160) 

del  (zet) 33  =  0  (161) 
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INPUT  VALUES  FOR  QUARTZ 

The  values  measured  by  Bechmann  (Ref.  1)  are  as  follows: 


TABLE 

6.  ISAGRIC 

ELASTIC 

STIFFNESSES. 

cEl  1 

cE12 

CE13 

CEl 4  CE33 

cE44 

cE66 

86.74 

6.98* 

11.91 

-17.91  107.2 

57.94 

39 . 88 

Units : 

10(9>  N/m2 

*  The  value  of  6.99  appearing  in  Ref.  1  has  been  changed  so 
that  the  relation  c66  =  (ell  -  cl2)/2  holds;  ell  and  c66  are 
directly  measured  and  hence  are  more  accurately  known  than 
cl2 . 


TABLE  7.  PIEZOELECTRIC  STRESS  CONSTANTS . 
ell  el4 

0.171  -0.0406 

Units:  C/m2 

TABLE  8.  DIELECTRIC  PERMITTIVITIES  AT  CONSTANT  STRAIN. 
(eps)Sll  (eps)S33 

39.21  41.03 

Units:  10(~12)  F/m 


OUTPUT  VALUES  FOR  QUARTZ 

The  input  values  from  Tables  6,  7,  and  8  were  used  to  compute  the 
remaining  elastic,  piezoelectric,  and  dielectric  quantities  for 
quartz  in  the  manner  discussed  in  prior  sections  of  this  report. 
The  results  are  given  in  Tables  9  to  16. 
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TABLE  9.  ELASTIC  STIFFNESSES. 


cE 

cD 

CP 

del  cDE 

del  cPE 

del  cPD 

11 

86.74 

87 . 49 

87.70 

0.746 

0.963 

0.218 

12 

6.98* 

6.23 

6.02 

-0.746 

-0.963 

-0.218 

13 

11.91 

11.91 

11.91 

0 

0 

0 

14 

-17 .91 

-18.09 

-18.14 

-0.177 

-0.229 

-0.0516 

33 

107 . 2 

107 . 2 

107 . 2 

0 

0 

0 

44 

57 . 94 

57.98 

57.99 

0.0420 

0.0543 

0.0123 

66 

Units : 

39.88 

10<9>  N/m2 

40.63 

40.84 

0.746 

0.963 

0.218 

*  The 

value  of  6. 

99  appearing  in 

Ref.  1  has  been  changed 

so 

that 

the  relation  c66  = 

(ell  - 

cl2)/2  holds; 

ell  and  c66  are 

directly  measured  and  hence  are  more  accurately  known  than 
cl2  . 


TABLE  10.  ELASTIC  COMPLIANCES. 


sE 

SD 

SP 

del  sED 

del  SEP 

del  sDP 

11 

12.77 

12 . 64 

12 . 60 

0.133 

0.171 

0.0379 

12 

-1.79 

-1.66 

-1.62 

-0. 133 

-0.171 

-0.0379 

13 

-1.22 

-1.22 

-1.22 

0 

0 

0 

14 

4 . 50 

4 .46 

4.45 

0.0419 

0 . 0538 

0.0119 

33 

9.60 

9 . 60 

9.60 

0 

0 

0 

44 

20.04 

20.03 

20.02 

0.0132 

0.0169 

0.00375 

66 

29. 12 

28.58 

28.43 

0.533 

0.684 

0. 152 

Units : 

10(-12)  m2/N 

TABLE  11. 

PIEZOELECTRIC 

[e],  [ h ] ,  AND  [a]  VALUES. 

e 

h  a 

11 

14 

0. 171 
-0.0406 

4.36  5.63 

-1.04  -1.34 

Units;  e: 

C/m2 ;  h  and  a: 

I 

I 

l 

I 

I 

l 

I 

l 

I 

I 

I 

l 

l 

l 

I 

I 

[ 

I 

1 

l 

e  i 

\  i 

>  i 
i 

O'  1 

O  1 

iH  | 

1 

1 
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TABLE  12. 

PIEZOELECTRIC  [d] ,  [g],  AND  [b] 

VALUES. 

d 

g 

b 

11 

2.31 

57.7 

74 . 1 

14 

+0.725 

18 . 1 

23.3 

Units:  d: 

io(‘12)  m/V;  g 

and  b:  10^-3^ 

m2/C 

TABLE  13. 

DIELECTRIC  (eps) 

VALUES . 

(eps)  S 

(eps) T 

del  (eps)TS 

11 

39.21 

39.97 

0 .759 

33 

41.03 

41.03 

0 

Units:  1 0 ( 

~12>  F/m. 

del  (eps)TS 

=  del  (chi)TS 

TABLE  14. 

DIELECTRIC  (chi) 

VALUES . 

(chi) S 

(chi)T 

del  (chi)TS 

11 

30.36 

31.12 

0.759 

33 

32 . 18 

32 . 18 

0 

Units:  10 ( 

~12'>  F/m. 

del  (chi)TS 

=  del  (eps)TS 

TABLE  15. 

DIELECTRIC  (bet) 

VALUES. 

(bet) S 

(bet)  T 

del  (bet) TS 

11  25.50  25.02  -0.485 

33  24.37  24.37  0 

Units:  10<9)  m/F. 
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TABLE  16.  DIELECTRIC  (zet)  VALUES. 


( zet) S 

( zet ) T 

del  (zet)TS 

11 

32 . 94 

32 . 14 

-0 .804 

33 

31.08 

31.08 

0 

Units : 

10(9)  m/F. 

CONCLUSIONS 

This  report  provides  formulas  interrelating  the  coefficients  that 
appear  in  the  several  alternative  sets  of  constitutive  equations 
involving  the  elastic,  piezoelectric,  and  dielectric  properties  of 
crystals.  These  are  then  specialized  for  crystals  of  class  32; 
using  measured  values  reported  for  quartz,  numerical  values  of  the 
elements  of  the  polarization  matrices  are  calculated. 
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White  Sands  Missile  Range,  NM 

221 

Cdr,  PM  JTFUSIOI! 

607 

Cdr,  Harry  Diamond  Laboratories 

ATTN:  JTF 

ATTN:  SLCHD-CO,  TO  (In  turn) 

1500  Planning  Research  Drive 

2800  Powder  Mi  1 1  Road 

001 

McLean.  VA  22102 

001 

Ade 1  phi  MO  20783-  1  145 

301  Rome  Air  Development  Center 

ATTN:  Documents  Library  (TILO) 
001  Griffiss  AFB,  NY  1344 L 


437  Deputy  for  Science  u  Technology 
Office,  Asst  Sec  Army  ( R&D ) 

001  Washington,  DC  20310 

438  HQDA  (OAMA-ARZ-O/Dr .  F.D.  Verderame) 

001  Washington,  DC  20310 

S20  Dir,  Electronic  Warfare/Reconnaissance 

Surveillance  and  Target  Acquisition  Ctr 
ATTN:  AMSEL-EW-0 

001  Fort  Monmouth,  NJ  07703-5000 

523  Oir,  Reconnaissance  Surveillance  and 
Target  Acquisition  Systems  Directorate 
ATTN:  AMSEL-EW-DR 

001  Fort  Monmouth,  NJ  07703-5000 

524  Cdr,  Marine  Corps  Liaison  Office 

ATTN:  AMSEL-LN-MC 

001  Fort  Monmouth,  NJ  07703-5000 

564  Oir,  US  Army  Signals  Warfare  Ctr 

ATTN:  AMSEL-SW-OS 
Vint  Hill  Farms  Station 
001  Warrenton,  VA  22186-5100 

602  Dir,  Night  Vision  1  Electro-Optics  Ctr 

CECGM 

ATTN:  AMSEL-NV-D 

001  Fort  8el voir ,  VA  22060-5677 


88002 


